Screening in (d + s)-wave superconductors: Application to Raman scattering 
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We study the polarization-dependent electronic Raman response of untwinned YBa2Cu307_£ superconduc- 
tors employing a tight-binding band structure with anisotropic hopping matrix parameters and a superconducting 
gap with a mixing of d- and s-wave symmetries. Using general arguments we find screening terms in the B lg 
scattering channel which are required by gauge invariance. As a result, we obtain a small but measurable soft- 
ening of the pair-breaking peak, whose position has been attributed for a long time to twice the superconducting 
gap maximum. Furthermore, we predict superconductivity-induced changes in the phonon line shapes that could 
provide a way to detect the isotropic s-wave admixture to the superconducting gap. 



I. INTRODUCTION 

The symmetry of the superconducting (SC) order param- 
eter in cuprate high-T c superconductors is now agreed to be 
unconventional after early intense theoretical and experimen- 
tal debates.[ 1] Phase-sensitive experiments such as corner- 
junction superconducting quantum interference device ex- 
periments [2] and tricrystal experiments on YBa2Cu307_a 
(YBCO) (Ref. H]) show that the SC order parameter under- 
goes sign changes in the Brillouin zone (BZ) that are con- 
sistent with the d x 2_ y 2 symmetry. A sign change of a SC 
order parameter with the d x i_ y i symmetry implies the ex- 
istence of four nodal points in the Brillouin zone that have 
been observed using momentum-resolved probes such as in- 
elastic neutron scattering (INS) or angle-resolved photoemis- 
sion spectroscopy (ARPES). Correspondingly, the simple gap 



A fc = AoCcosfea, -COSfcj,)/2 



(1.1) 



with a d x 2_ y 2 symmetry on the Brillouin zone has often been 
used as a starting point of a quantitative interpretation of these 
experiments. 

Of course, a sign change in the SC gap does not preclude a 
gap that is more complicated than the simple d x 2_ y 2 gap ( 11.11 ). 
For those high-T c cuprates with a tetragonal crystalline struc- 
ture, higher rf-wave harmonics have been invoked to explain 
ARPES measurements of the magnitude of the SC gap.|0] 
For the cuprate family YBa2Qi307_,5, which exhibits quite 
strong orthorhombic distortions, one expects corrections to 
the d x 2_ y 2 gap ( 11.11 ) on symmetry grounds alone. This expec- 
tation has been confirmed by several experimental methods 
such as ARPES studies,J5|]INS studies, (6D and measurements 
of Josephson currents, |7j l8[] and has also been investigated 
theoretically.HHEl] 

Polarization-dependent electronic Raman scattering also 
probes the momentum dependence of the magnitude of the su- 
perconducting order parameter and has provided yet one more 
piece of evidence for the d x 2_ y 2-wave pairing scenario. 
[l3l[3 In particular, for tetragonal high-T c cuprates, one finds 
(a) various low-energy power laws in different polarization 
channels that are consistent with the existence of nodal points 



for the gap and (b) the pair-breaking peak in the B lg channel 

at energy twice the superconducting gap maximum A seen 
by other means. fl5ll 

In this paper, we are going to investigate the consequences 
for polarization-dependent Raman scattering of a subdom- 
inant admixture of an isotropic s-wave component to the 
gap ( 11.11 ) which should be of relevance to orthorhombic 
high-r c superconductors of the YBCO family. Assuming 
the existence of well-defined SC quasiparticles, we compute 
the polarization-dependent electronic Raman-scattering cross 
section including the effects of (i) an orthorhombic tight- 
binding dispersion with (d+s) gap, (ii) a long-range Coulomb 
interaction treated within the random-phase approximation 
(RPA), (iii) and the effective mass approximation for the Ra- 
man vertex. We show that the pair-breaking peak in the B lg 
channel is softened by an amount proportional to the isotropic 
s component to the gap in that it occurs at a lower value than 
the absolute maximum of the gap on the "normal-state" Fermi 
surface. We also show that the A lg channel develops a dou- 
ble peak structure with the peak separation proportional to the 
isotropic s component to the gap. Furthermore, we compute 
superconductivity-induced changes in the phonon line shapes, 
and argue how Raman scattering on phonons allows us to ex- 
tract a signature of a subdominant and isotropic s-wave com- 
ponent to the gap. 



II. THEORY 

The differential cross section in a Raman-scattering exper- 
iment for a momentum transfer q that is small compared to 
the extension of the Brillouin zone is proportional to [1 + 
n(w)]x"(w), where n denotes the Bose distribution, u> the 
frequency of the incoming plane wave, and 

X 7 H = ( x ; + ix") M = X 7 (9 « 0, u + irf) (2.1a) 
is the linear-response function for the density operator 



EE 

fe <r=T4 



(2.1b) 
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The coupling between the SC quasiparticles, linear superpo- 
sitions of the fermionic creation ct and annihilation c k a 
operators, and the incoming (outgoing) plane wave with the 
polarization vector e 1 (e ) is here approximated by 

^E^J^ (2.10 

in the nonresonant limit. The "normal-state" dispersion 

£fc = -2* [(l + S )cosk x + (l-6 )cosk y ] 

—At' cos fc x cos k y — (i (2. Id) 

(see Figs. 1 and 2 in Ref. Ilq) combines with the gap 
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A fe = -y (cos fc^ - cos fej,) + A s 



(2.1e) FIG. 1: Calculated electronic Raman response of (d + s)-wave su- 
perconductors for the -Big, B 2g , and A lg channels (tetragonal no- 



to give the SC quasiparticle dispersion E k = e| 
Both <5 and A s represent symmetry -breaking terms that lower 
the symmetry from tetragonal to orthorhombic in an effective 
one-band description of a single copper-oxygen plane. The s- 
wave component A s is isotropic [compared with the extended 
s-wave admixture from Eq. (13.81) 1. The Raman vertices j k 
can be classified according to the irreducible representations 
of the symmetry group of the crystal. For a crystal with tetrag- 
onal symmetry (point group D^A, the relevant symmetries are 
the B lg , B 2g , and A lg polarizations. As we shall consider a 
model with subdominant orthorhombic distortions, 5 <C 1, 
A s <C A , we will use, in what follows, the notation of tetrag- 
onal symmetry. 1 1711 If so, we can identify the B 1 , B 2g , and 
A x channels for Raman scattering with the Raman vertices 



tation), respectively, according to Eq. q273j. The B lg pair-breaking 
peak softens to u> ~ 44meV < 2A ~ 56meV where A is the SC 
gap maximum on the FS denned in Eq. P. It . The A lg spectrum has 
been multiplied by the factor of 20. Inset: "bare" density of states in 
the superconducting state for orthorhombic symmetry. 



in the B lg , B 2g , and A lg channels for Raman scattering, re- 
spectively. As usual, the bracket (in a box of volume V) 



((•••)<U, 



k 



• )A| tanh (ft) (2.4) 



VIE,. 



in-2E h 



denotes an average over the BZ weighted by the Tsuneto func- 



lB lg fc oc t [(1 + 5 ) cos k x -{l- S ) cos k y ] , (2.2a) tion k .[22] The second term in Eqs. (|2.3a|) and (j2.3c|) is com 

7b 2 k x 4t' sin k x sin k y , (2.2b) 
lA lg k « t [(I + S ) cos k x + (I -5 ) COS k y ] 

+At' cos k x cos k , (2.2c) 



respectively. 

The electronic Raman response is calculated in the gauge 
invariant form assuming that the quasiparticles interact 
through the long-range Coulomb potential V q = ^p- . For the 
tetragonal symmetry, the RPA for the polarization-dependent 
Raman resp onse function x 7 ( w ) was derived in Refs. fl~5l fTsI 
and [14J, lljj |2fj|, |2lj and shown to be well defined in the limit 

Its generalization to orthorhombic symmetry leads to 
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-TbJh 



Xa 1s M = (lljk) 



lA^k 



(2.3a) 
(2.3b) 

(2.3c) 



monly called screening term. It can be viewed as originating 
from the Goldstone mode of the superconductor and ensures 
gauge invariance of the Raman response. II 1511 For tetragonal 
symmetry only screening terms in the A lg scattering channel 
are possible. That is, the screening term in Eq. d2.3ab is van- 
ishing identically if the SC quasiparticle dispersion is of the 
tetragonal symmetry. In the presence of orthorhombic distor- 
tions of the YBCO type screening terms can affect the B lg 
channels on general symmetry grounds, but are absent in the 
B 2g channel. Similarly, in orthorhombic Bi2Sr2CaCu208+£c, 
whose crystallographic axes are rotated by 45° with respect 
to the CuO bonds, screening terms arise in the B 2g channel, 
but do not affect the B lg channel. In the following we shall 
only consider orthorhombic superconductors of the YBCO 
type and disregard any possible screening terms in the B 2g 
channel. The presence of screening terms in the B lg or B 2g 
channels for orthorhombic superconductors has been previ- 
ously reported in the literature, see Refs.lldandl2lL 



III. DISCUSSION 

We have computed numerically the Raman response d2.3l ) 
using the parameters t = 200 meV, 5 a = -0.03, t'/t = -0.4, 
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lift = -1.2, A = 30 meV, and A s = 6 meV. We 
also choose the damping 77 = 1 meV. This choice yields a 
SC quasiparticle dispersion that is in qualitative agreement 
with the one measured by photoemission experiments on un- 
twinned YBCO close to optimal doping. j23l 

In Fig.Q]we show the resulting Raman response for (d + s)- 
wave superconductors. We find that the low-energy asymp- 

totics, Xs is M oc (w/2A„) , Xs 2s M « (w/2A ), and 

x'a ( w ) k ( W /2A ^ , remain essentially unchanged with or 
without the symmetry -breaking terms S and A s . Here we 
have introduced the absolute maximum 



\ n = max A(fe) 
fceFS 



28meV 



(3.1) 



of the gap over the "normal-state" Fermi surface (FS) with full 
tetragonal symmetry (<5 = A s = 0). We also find that the rel- 
ative peak heights in the B lg and A lg responses remain lfl5ll 
after switching on S Q > and A s > 0. These properties are 
robust to the lowering of the tetragonal to the orthorhombic 
symmetry as they owe their existence to that of the nodes of 
the SC gap. There are nevertheless important changes in the 
line shapes of the A lg and B lg responses due to orthorhom- 
bicity which reflect the change in the underlying density of 
states (see inset of Fig. Q]). Most importantly, the B lg pair- 
breaking peak is not located any longer at the absolute gap 
maximum lu = 2A ~ 56 meV. Instead, it is shifted by 
~ 2A S to a lower value due to the additional screening term in 
Eq. ( I2.3ab and thus its position tracks the position of the dou- 
ble peaks in the A lg channel. B24I1 This is interesting because it 
demonstrates that the common lore that the B lg pair-breaking 

peak is always at w = 2A is not correct in (d + s)-wave su- 
perconductors. In passing, we note that a comparison between 
the (XX)-polarization channel, e° = e 1 = (J), and (re- 
polarization channel, e° = e 1 = (?), reveals only minor 
differences in the peak positions and line shapes (not shown). 

Recently, it has been proposed that vertex corrections to 
the Raman response function due to a spin-mediated pairing 
interaction can lead to the development of a resonance be- 
low 2A n 



2a 1221 |28fl Inclusion of these interaction effects 
would lead to an additional shift of the B 1 pair-breaking peak 
to lower frequencies. However, the magnitude of this shift de- 
pends on model assumptions and it remains an open question, 
whether these vertex corrections are of any relevance for the 
interpretation of the Raman data on high-T c cuprates. 

In order to analyze in more detail the softening of the B lg 
pair-breaking peak due to the screening term, we show in 
Fig. [2] the results for A s = meV, 3 meV, and 6 meV, re- 
spectively, with all other parameters held fixed. We find that 
increasing A s induces a reduction of the peak intensity and 
a shift to lower frequencies of the peak positions. A closer 
inspection of Eq. d2.3a| ) reveals that the bare contribution 



B 



(3.2) 



exhibits two peaks located at u> = 2(A =p A s ), respectively. 
This double peak structure is a reflection of the bare quasipar- 
ticle density of states that features two local maxima (inset of 
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FIG. 2: (color online) Calculated electronic Raman response of [d + 
s)-wave superconductors for the B lg channel (tetragonal notation) 
according to Eq. J2.3b . The fermiology parameters from Eq. d2.lt 
are the same as in Fig.fTJexcept for A s = meV (red), A s = 3 meV 
(blue), and A s — 6 meV (black). Inset: decomposition of the left- 
hand side of Eq. d2,3a| > into the sum of the bare contribution Bi g and 
the screening contribution Bf^ on the right-hand side of Eq. J2.3at 
when A, = 6 meV. 



Fig-IU) at half the frequencies of the peaks in B® g . In contrast, 
the screening contribution 



B^{uj) = Im 



7b, OS /(9 k ) u 



(3.3) 



shows a broad dip followed by a peak (see the inset of Fig.f2|. 
That is, an enhancement of the scattering efficiency by screen- 
ing {antiscreening) occurs at low frequencies u> < 2Ao, 
whereas at high frequencies w > 2 Ao the scattering efficiency 
is reduced by the screening term. To show how antiscreening 
arises, we rewrite the imaginary part of the B lg screening term 



nscr 



M = 



7b, 



(3.4) 



7s, 



7b, 



((>x + to 

where the real and imaginary parts of the Tsuneto functions 
are denoted by 0' k and 6'^, respectively. It is instructive to 
approximate the quasiparticle dispersion by 



E, 



k,<j> 



|A„ 



(3.5a) 



with the linearized dispersion e k = v F (k — k F ) and the gap 
function A^ = A cos(2</>) + A s that depends only on the 
Fermi surface angle <f>. Similarly, we replace the B lg vertex 
by its angular-dependent form, 



7b, 



cos(20). 



(3.5b) 
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With these simplifications and in the zero-temperature limit 
the imaginary part of the BZ averages reduce to 



where we have restricted the BZ integration to a ring — 8k < 
k — k F < 5k around the Fermi surface, the summation runs 



over all <A, £ 



lo = 2E 



h,(j> 



and A'i denotes the deriva- 



tive of the gap function A^. The imaginary part of the av- 
eraged Tsuneto function {9'D is a positive function of lo and 
exhibits a positive divergence whenever A'^ = 0, i.e., at the 
frequencies 



„( 2 ) 



2(A -AJ, 



E, n = 2(A + AJ, 



(3.7) 



as can be seen from the last line of Eq. i3.6i . On the other 
hand, (7^ 6") possesses a negative divergence at a4 and a 

positive divergence at Wc , and changes its sign at ~ 2A , 
since the B lg vertex j B ~ cos 2<f> exhibits a sign change 
along the Fermi surface. Any divergence in the frequency de- 
pendence of (9") and (7^ 9") results in a steplike disconti- 

nuity in (9 1 ) and (7^ 9'), respectively, due to the Kramers- 
Kronig relation. Vice versa, from the absence of any step- 
like discontinuity in the imaginary part of the BZ averages 
follows the absence of any divergence in the real part of 
the corresponding quantity. Furthermore, it turns out that 
both {9') and (~f B 9') are positive in the frequency range 

< lo < 2(Ao + A s ). Taking all these observations together, 
we find that the first term in the B lg screening function d3.4b 

features a negative divergence at lo = 2(Ao — A s ), which 
is compensated by a positive divergence of the second term. 
But then, at lo = 2(Ao + A s ) both terms in Eq. (13.41 i show a 
positive divergence, which results in a strong screening of the 
second peak in B® . At lo ~ 2Ao the first term in Eq. d3.4| > is 
vanishing, whereas the second term is negative, which yields 
to antiscreening of the B lg scattering efficiency. 

It is important to note that the effects induced by a nonvan- 
ishing isotropic s-wave gap A s in Fig.[2]would not occur had 
we assumed a subdominant extended s-wave gap such as 



A fc = A o (coak x 



cos k y ) /2 



A° xt (cos fes + cos/y /2. 



(3.8) 



Contrary to the isotropic subdominant s-wave admixture A s 
in Eq. ( 12. let , the extended s-wave admixture in Eq. ( 13.8b only 
leads to a shift of the nodal line, but does not give a different 
absolute gap value at (it, 0) compared to (0, ir) (see Fig. 2(d) 
in Ref. 16). From the above discussion one infers that the 
property |A^ =? | ^ | A^ =0 | is crucial for the splitting of the 



B lg pair-breaking peak as well as for the appearance of antis- 
creening in the B lg channel. Thus, we conclude that to a first 
approximation an extended s-wave admixture leaves the B lg 
response unchanged, an expectation that we have confirmed 
by numerical calculations. 



IV. SUPERCONDUCTIVITY-INDUCED CHANGES IN 
PHONON LINE SHAPES 

In principle, the existence of a subdominant and isotropic 
s-wave component A s to the rf-wave SC gap can be extracted 
from a line-shape analysis of the A lg and B lg electronic re- 
sponses in a polarization-resolved Raman-scattering experi- 
ment. However, as part of the electronic Raman signal is in 
general masked by phonon excitations, it might be hard to 
detect experimentally these changes in the line shape. On 
the other hand, we argue in the following that polarization- 
resolved Raman scattering on phonons can be used to extract 
a signature of A s . 

To substantiate this point, we calculate the 
superconductivity-induced changes in the self-energy of 
optical, zone-center (q — 0) phonons. Thereto, we consider a 
linear coupling between electrons and phonons 



cl— ph 



1 

v 



E 9l 



q C k+q,a C k,a 



1,1 



(4.1) 

where b q and b q are the creation and annihilation opera- 
tors of phonons in a given branch 7 with phonon frequency 
y , respectively, and denotes the electron-phonon cou 



pling. The form of the electron-phonon interaction d4.il ) is 
the one for the nonresonant electronic Raman scattering pro- 
vided the effective Raman vertex is replaced by the electron- 
phonon vertex . Hence, within an RPA treatment of the 
Coulomb interactions, we find that in the q — > limit the 
superconductivity-induced changes in the phonon self-energy 
are given by iiM IH 



s» = -(UZ.o 



(4.2) 



where the angular brackets are defined by Eq. ( 12.4b - The sym- 
metry of the optical phonons is encoded in the matrix ele- 
ment gjl . The electron-phonon coupling of the A lg and B lg 
phonons are in a first approximation given by 



9ko=9b, (cosk x - cosfc )/2, 



9hfi = 9A lg (cos k x + cos k y )/2. 



(4.3) 



with the electron-phonon coupling constants g B and 
g A . In Fig. [3] we have numerically evaluated the 

superconductivity-induced changes in the self-energy S 7 (w) 
for B lg and A lg phonons in a (d + s)-wave superconductor. 
To estimate the size of these effects we have assumed some 
typical values for the phonon coupling strength g 7 that are in 
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FIG. 3: (color online) Panel (a) displays the real part of the 
superconductivity-induced phonon self-energy ReE Bi (o>) of the 

B lg phonon with g% i = 0.1A /Nf and at zero temperature ac- 
cording to Eq. i4.2l . The fermiology parameters are the same as 
in Fig. 1 except for A s , which varies form A s = meV (red), to 
A s = 3 meV (blue), and A a = 6 meV (black). The inset shows 
the decomposition of the left-hand side of Eq. \4.2l into a sum of the 
bare contribution ReEo and the screening contribution ReEg r 

1.9 ig 

when A s = 6 meV. Panel (b) shows the real and imaginary parts of 
the superconductivity-induced self-energy of the A lg phonon with 
g% = 1.0A /Nf and the same fermiology parameters as in panel 

(a) . Negative (positive) values correspond to softening (hardening) 
in the main frame and broadening (sharpening) in the inset of panel 

(b) . Here Nf denotes the "normal-state" density of states per spin at 
the Fermi energy. 



rough overall a gree ment with the observed Raman shifts in 
YBCO.iH Ull iiA negative (positive) Re S 7 (w) 
corresponds to a softening (hardening) of the phonon below 
T c , whereas a negative (positive) ImS 7 (w) leads to a broad- 
ening (sharpening) of the phonon linewidth below T c . 

For B lg phonons [see Fig. [3fa)] we find that the real part 

of the phonon self-energy crosses zero around 2Ao irrespec- 
tive of the value of A s . However, upon inclusion of a sub- 
dominant s-wave gap the maxima and minima of RcE R are 

is 

shifted to 2(Ao + A s ) and 2(Ao — A s ), respectively. Hence, 
the maximal softening of the phonon frequency occurs at a 



frequency 2A S smaller than the gap maximum Ao. The in- 
set of Fig. |3 a) displays the decomposition of the real part of 
the self-energy Re E B (lo) into its screened Re E^ 1 (10) and 

bare parts ReEn (uo). The frequency dependence of these 

two contributions mimics that of B® g (oj) and Bl g r (oj), which 
we discussed in the previous Section. A finite s-wave com- 
ponent splits the step in Re £ ° B (uj) into two steps located at 

2(Ao =F A s ) . In between these two steps the screening con- 
tribution Bl c g r (ui) is negative, which leads to antiscreening. 
Neglecting the influence of a nonzero <5o, one can infer the 
Kramers-Kronig transform of the curves in Fig. [2a), directly 
from Fig. [2] and thereby obtain the imaginary part of the self- 
energy ImS B (uj). Due to the s-wave admixture, the fre- 
quency where the maximal broadening occurs is shifted from 
2A to 2(A - AJ.JH 

Similar effects occur for the self-energy of A lg phonons 
[see Fig.(3fb)]. An isotropic s-wave component splits the step 
in the real part Re (oj) into two steps, thereby shifting 

the crossing point of Re T^ A ^ (uj) with the zero line to lower 
frequencies (A s = 3 meV) or to higher frequencies (A s = 
6 meV). Contrary to B lg phonons, antiscreening effects are 
absent in the self-energy of A lg phonons. 

Finally, we note that the effects on the superconductivity- 
induced changes of the phonon self-energy induced by an 
isotropic s-wave gap are again absent for a gap of pure d- 
wave character or for an extended s-wave admixture such as 
in Eq. ( 13.81 ). The effects on the phonon self-energy reported 
in this paper therefore constitute a fingerprint of an isotropic 
s-wave admixture to the pairing symmetry. 



V. CONCLUSIONS 

We have calculated the polarization-dependent Raman re- 
sponse for electrons and phonons in (d + s)-wave supercon- 
ductors. We find that the screening terms in the B lg channel 
lead to a softening of the B lg pair-breaking peak of the order 
of 2A S , i.e., twice the value of the isotropic s-wave compo- 
nent to the SC gap. This fact calls into question the long- 
standing interpretation that the B lg pair-breaking peak is lo- 
cated at oj = 2A and thus can be directly used as a measure 
of the rf-wave component of the SC gap when orthorhombicity 
is present. Secondly, we have estimated the effects of a sub- 
dominant s-wave admixture on the superconductivity-induced 
phonon renormalizations. These effects, although compara- 
tively small, might serve as a fingerprint of an isotropic s- 
wave admixture to the pairing symmetry. 
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